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ABSTRACT 

We report the recovery of a spectroscopic event in rj Carinae in 1997/98 
after a prediction by Damineli (1996). A true periodicity with P = 2020 ± 5 
days (0.2% uncertainty) is obtained. The line intensities and the radial-velocity 
curve display a phase-locked behavior implying that the energy and dynamics 
of the event repeat from cycle to cycle. This rules out S Doradus oscillation or 
multiple shell ejection by an unstable star as the explanation of the spectroscopic 
events. A colliding-wind binary scenario is supported by our spectroscopic data 
and by X-ray observations. Although deviations from a simple case exist around 
periastron, intensive monitoring during the next event (mid 2003) will be crucial 
to the understanding of the system. 

Subject headings: stars: individual (rj Carinae) - stars: variables: other 
(luminous blue variable) - stars: emission line - stars: binary 

1. Introduction 

Eta Carinae is a key object to understand the evolution of massive stars. It is 
one of the most luminous stars known and underwent a great eruption that created the 
Homunculus 150 years ago. Since then, the nature of rj Carinae has been the subject of an 
intense debate. However, the star does not fit any of the proposed models - see Davidson 
& Humphreys (1997) for a review. In the last three decades the prevailing view was that rj 
Car is a single luminous blue variable star (LBV). The spectroscopic events - fading of high 
excitation lines - occasionally observed, were interpreted as being imprints of S Doradus 
oscillations (Damineli et al. 1998). The regularity of the 5.5 year cycle, however, makes 
such a hypothesis unlikely (Damineli 1996). The binary system proposed by Damineli, 
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Conti & Lopes (1997 - hereafter DCL) brought new insight to the problem. The solution 
to the radial-velocity (RV) curve of the broad component in the Paschen— 7 line resulted in 
two massive stars {M 1 ~M 2 ~ 70M Q ). The large mass-loss rates of the stars would imply 
strong wind-wind collision producing thermal hard X-rays. The highly eccentric orbit would 
produce phase-locked X-ray variability. Those features are in accord with the observations 
(Corcoran et al. 1998) indicating that 77 Car is a colliding-wind binary (CWB). Davidson 
(1997) presented an alternative solution to the DCL data, leading to similar characteristics 
but larger eccentricity. The parameters of the binary system, however, were not trustworthy 
because they were derived from only one cycle and based on a single emission line. The 
prediction of a spectroscopic event for late 1997 through early 1998 opened an outstanding 
opportunity to test this model against the, at the time dominant, idea of recurrent shell 
ejection. 

The expected event indeed occurred at the right time as shown by the spec- 
troscopic data analyzed in the next section. The event was also detected at radio 
wavelengths (Abraham & Damineli 1999) and in X-rays, displaying a dramatic variation 
(Ishibashi et al. 1999). A full 2-D hydrodynamic simulation (Stevens & Pittard 1999) 
and even a simple CWB model (Ishibashi et al. 1999) reproduce the overall behavior of 
i] Car in X-rays. However, significant discrepancies remain around periastron. This is 
not unexpected, taking into account that near periastron the stars are close enough for 
deviations from spherical symmetry to be important. The wind is probably non-spherical 
and an equatorial disk may surround the primary star. 

The event observed in optical/near- infrared spectral lines and at radio wavelengths is 
due to a cause different from that in X-rays, despite their temporal coincidence. Following 
Damineli, Lopes, & Conti (1999), the secondary star is the main source of hard photons, not 
the wind-wind colliding zone. The immersion of the secondary star into the companion's 
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wind prevents hard photons from reaching the external regions of the wind of the primary 
and the circumstellar gas. The resulting effect is a temporary drop in the degree of 
ionization and, consequently, the fading of the high excitation lines. 

Although preliminary results indicate that rj Car is likely to be a binary, critical tests 
must be carried out to provide a definitive proof. Such tests are provided in this letter, 
based on the data collected in a spectroscopic campaign initiated in 1989 that covered twice 
the 5.5 year cycle. 

2. A True Periodicity 

Spectroscopic observations of rj Car were carried out at the European Southern 
Observatory - ESO/Chile and at the National Astrophysical Laboratory - LNA/Brazil 
- see Wolf et al. (1999) and Damineli et al. (1999) for details. The [S ill] and 
[Anil] spectral lines plummeted in November 1997 and disappeared a month later 
(Lopes & Damineli 1997). Around that time, the near-infrared light-curve reached a local 
maximum (Whitelock & Laney 1999) and the spectroscopic lines a minimum, similar to 
those reported by Whitelock et al.(1994) and by Damineli (1996). Since we have detailed 
observations of the 1992 event, we folded them with the 1998 event. The sawtooth structure 
in the RV curve provides an accurate measure of the elapsed time between the last two 
events. We plotted in Figure la the RV variations of the broad component in He I A6678 for 
the 1998 event (solid line) superposed with the previous event shifted by 2020 days (dotted 
line). The matching was obtained by shifting one relatively to the other by 2020 ± 5 days. 
The same recurrence time was derived from Pa-7 although the uncertainty is larger for this 
line due to a smaller number of observations. 

We should not expect, a priori, that line intensities would result in the same elapsed 
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time as for RV. The narrow component in [S in] A6312 and He I A6678 lines disappeared 
around 12 December 1997 (JD 2450794±2 days). The same was observed on 2 June 1992 
(JD 2448775±5 days), resulting in a recurrence time 2019 ± 7 days, compatible with that 
from RV curve. This time scale is also in agreement with P = 2014 ± 50 days proposed by 
DCL based on the 1948, 1965, 1981, 1987 and 1992 events. This is an indication that the 
previous events resembled closely the last two ones we have observed: deep and brief. If the 
events had lasted longer, no obvious strict periodicity would have been hinted from such a 
short list of occurrences. If line intensity fadings had not been so remarkable they would 
have been mistaken as secondary fluctuations that are frequent in the rj Car spectrum. 

An additional way of checking whether the periodicity is true is by comparing our 
data with those of Gaviola (1953). The description of the spectrum collected by that 
author in 1948 is detailed enough to show that rj Car was very close to the center of a 
spectroscopic event on April 19th. He reported that [Nil] A5755 was fainter than [Fell] 
A5748. Our spectra collected during the last two events show that such an inversion of 
intensities remains for less than three months around the center of an event. The 9 cycles 
between the 1948 and the 1998 events constrain very well the uncertainty in the periodicity: 
P = 2019 ± 10 days. This figure is in excellent agreement with that described above, based 
on the last two events. We can conclude that the events in rj Car are truly periodic with 
P = 2020 ± 5 days (or P = 5.53 ± 0.01 years). 

The equivalent widths of the spectral lines, especially He I and the high-excitation 
forbidden lines, are quite predictable from cycle to cycle. In Figure lb we display the 
He I A6678 line during the events of 1998 and 1992 (the latter one shifted by 2020 days). 
This figure implies that the emitting volume and the speed of gas display a phase-locked 
behavior. Due to the fact that no luminous hot star has been observed pulsating in such a 
regular fashion and ejecting identical shells, we rule out stellar instability as the mechanism 
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responsible for the spectroscopic events in r/ Carinae. 

3. The Orbital Parameters 

The orbital solution presented by DCL had been discussed by Davidson (1997, 1999). 
Some concerns have been raised about the reliability of the broad-line components to trace 
the orbital motion. On the one hand, broad lines are formed inside the wind and should 
display Doppler shifts if the star were a binary. On the other hand, radiative-transfer effects 
could mask Doppler motion. If the event is produced by the plunging of the secondary star 
into the primary-star wind (Damineli, Lopes, & Conti 1999), the ionization structure of 
the wind should change around periastron, mixing together Doppler and radiative-transfer 
effects. Different lines would display different line-profile variations, making it difficult to 
predict what line tracks the orbital motion better. The best candidates are lines (broad 
components) displaying small changes in equivalent width and in line shape, with as faint 
as possible P Cygni absorption components. Some hydrogen Paschen lines are free of blends 
and look suitable for this purpose, as seen in DCL (their Figure 1). 

He I lines show strong variability in intensity and line-profile. Deblending the broad 
components from the nebular and P Cygni contamination is not as straightforward as in 
H-Paschen lines. The observed changes in He I lines reflect more the excitation effects in 
the wind than Doppler motion of the stars and this is why we did not use them in the RV 
solution. The large number of observations we collected for He I A6678 line, however, still 
makes it the best tool to measure the period length. 

The RV curve containing all the available Pa-7 and Pa-5 data results in the solution 
displayed in Figure 2. The best fit gives: e = 0.75, u>i = 275°, T periastron =JD2450861 (17 
February 1998), K 1 = 50 km s -1 and /mi = 7.5M . Conjunctions occur at phases 0.999 
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and 0.438. The standard deviation is a — 11 km s _1 . This result is in general agreement 
with those of Damineli, Lopes, & Conti (1999) and Davidson (1997). The main difference 
between present and previous results is the smaller mass function, (/mi), accommodating 
a smaller mass for the secondary star, which is desirable, since its spectrum is not visible in 
the optical or the UV. The mass of the unseen companion is difficult to be constrained, as 
it depends on the mass of the primary-star, the mass-loss rate and the orbital inclination. 
Regarding the primary-star we adopt DCL's Mi ~ 70M Q , based on the total luminosity of 
the star, the age of Trl6 cluster and standard evolutionary models. 

Figure 2 gives strong support in favor of binarity in rj Car, but it does not necessarily 
guarantee that our orbital elements are accurate. There are two problems with our RV 
curve. First, Pa-7 shows changes in the shape of the broad-line component, indicating that 
the wind is not rigidly following the star. The effect in the RV curve seems to be small 
because the 7 velocity (-12 km s~ x ) derived from the RV curve agrees with V mrf = — 7 ± 9 
km s -1 obtained by Conti, Leep, & Lorre (1977) from seven O-type stars in the Carina 
Nebula. Second, data around the periastron are scarce, which is a problem for a system 
with such a high eccentricity. Acceptable eccentricities range from e~0.65 up to e~0.85, 
displacing the time of periastron by several months (up to 0.05 in phase) around the center 
of the spectroscopic event. The other orbital elements are less sensitive to the details of the 
RV curve around periastron. 

4. Discussion 

The recovery of the 1998 spectroscopic event within <1% of the predicted time 
supports the hypothesis of true periodicity. The strict repeatability in the RV curve and 
line intensities is unambiguously in favor of the binary model and against the idea of an 
unstable star as the explanation of the 5.5 year cycle. Regardless of the particular value 
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adopted for the orbital eccentricity (e>0.65), the primary star is a typical LBV and the 
unseen companion a hotter and less evolved star. The orbital parameters must be regarded 
as very provisional, because data around periastron are scarce and based on broad emission 
lines and the physical parameters of very massive stars are uncertain. Nevertheless, the 
success of the CWB models in reproducing X-ray luminosity, temperature, Nh, and the dip 
in the light curve, add credibility to the solution. Fitting the X-ray light curve by CWB 
models presents some problems around periastron, as shown by Ishibashi et al. (1999) 
and Stevens & Pittard (1999). However, non-spherical symmetry in the stellar wind or a 
circumstellar disk might remove the disagreement between observations and CWB models. 

However, the rj Car wind is not steady and probably not homogeneous, producing a 
RV curve not as simple as in classical cases. Significant profile variations have been seen 
in some spectral lines around periastron (Davidson 1999), indicating that gas blobs are 
thrown away from the stars, a reminiscence of shell ejection suggested in previous works. 
This time, however, the mechanism is not advocated to explain the whole 5.5 year cycle 
behavior, but minor perturbations. 

Although the binary scenario accounts for the main observational characteristics of 
7] Carinae, the nature of the spectroscopic events is not fully understood. Continuous 
monitoring of the system, with a better time sampling than the previously attained, is 
needed for the next event in 2003. An observational campaign starting in May should 
extend until late September to determine the RV curve, to late October to survey the dip in 
the X-ray light-curve, and, through December to follow the variability in different spectral 
lines. The star is not accessible during night-time from August through October, which 
will preclude ground-based optical observations. High spatial resolution will be crucial to 
disentangle the stellar intrinsic variability from the associated light reflections through the 
surrounding nebula. Far-UV spectra are relevant, as the secondary (hotter) component of 
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the system may be detectable; consequently, space-based facilities will be invaluable. From 
the ground, infrared and near infrared will be required to perform daylight observations. 

We thank the ESO and LNA Observatories for the generous allocation of observing time. 
This work was supported by FAPESP, FINEP and by the Deutsche Forschungsgemeinschaft 
(Wo 296/20-1). A. D. thanks G. Medina Tamco for helpful discussions. 
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Fig. 1. — Variability during the present cycle (solid line) folded with the previous one 
displaced by 2020 days (dotted line); a) radial velocity (RV) of the broad component in the 
Hel A6678 emission line; b) equivalent width (EqW) of the blend He I A6678+[Nill] A6666. 
Dotted tickmarks indicate beginning of years. 
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Fig. 2. — Radial velocity curve of Pa— 7 and Pa— 5 broad-line components. Previous cycle 
(circles), present cycle (crosses), and orbital solution (solid line). 



